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ABSTRACT 



We present spectroscopic observations of a sample of 15 embedded young 
stellar objects (YSOs) in the Large Magellanic Cloud (LMC). These observa- 
tions were obtained with the Spitzer Infrared Spectrograph (IRS) as part of the 
SAGE-Spec Legacy program. We analyze the two prominent ice bands in the IRS 
spectral range: the bending mode of CO2 ice at 15.2 fj,m and the ice band between 
5 and 7 /im that includes contributions from the bending mode of water ice at 
6//m amongst other ice species. The 5— 7//m band is difficult to identify in our 
LMC sample due to the conspicuous presence of PAH emission superimposed 
onto the ice spectra. We identify water ice in the spectra of two sources; the 
spectrum of one of those sources also exhibits the 6.8-//m ice feature attributed 
in the literature to ammonium and methanol. We model the CO2 band in de- 
tail, using the combination of laboratory ice profiles available in the literature. 
We find that a significant fraction (> 50%) of CO2 ice is locked in a water-rich 
component, consistent with what is observed for Galactic sources. The majority 
of the sources in the LMC also require a pure-C02 contribution to the ice pro- 
file, evidence of thermal processing. There is a suggestion that CO2 production 
might be enhanced in the LMC, but the size of the available sample precludes 
firmer conclusions. We place our results in the context of the star formation 
environment in the LMC. 

Subject headings: astrochemistry — circumstellar matter — stars: formation — 
galaxies: individual (LMC) — Magellanic Clouds 



- 3 - 



Introduction 



The formation of stars in the high-redshift Universe occurred in a metal-poor 
environment. Much of what we know about the process of star formation is deduced from 
observations in nearby, Galactic star-forming regions with essentially solar metallicity. 
However, many of the physical processes that affect star formation and their observational 
diagnostics are expected to depend on metallicity, for instance cooling and gas-phase and 
grain surface chemistry. The nearest templates for the detailed study of star formation 
under metal-poor conditions are the Magellanic Clouds (MCs). The Interstellar Media 
(ISM) of the Large and Small Magellanic Cloud (LMC and SMC) have significantly lower 
metalhcities (.^lmc ~ 0.4 Z0 and Zsmc ~ 0.2 Z©), than the solar neighborhood ISM {Zcai ~ 
Zq). Only relatively recently are we able to s tudy the deta ils of the star formation process 



in the MCs (for a review of recent results see 



Oliveira 



2OO9I). 



The star-formation process is a complex interplay of various chemo-physical processes. 
During the onset of gravitational collapse of a molecular cloud, sufficiently dense cores 
can only develop if the heat produced during the contraction can be dissipated. The most 
efficient cooling mechanisms are via radiation through fine-structure lines of carbon and 
oxygen (e.g., via the emission lines of [Cll] at 158 /xm and [Ol] at 63 and 146 um), and 



rotat ional transitions in abundant molecules such as CO, O2 and water (IGoldsmith fc Langei 
IQTsI ). These cooling agents all contain at least one heavy atom — H2 has no permanent 
dipole moment, and emission through its rotational transitions is very unlikely, making it 
an inefficient coolant. Besides also contributing to the thermal balance in the molecular 
cloud, dust grains are crucial in driving cloud chemistry, as dust opacity shields molecules 
from radiation and grain surfaces enable chemical reactions to occur that would not 
happen in the gas phase. Surface chemis try is also thou ght to play an importan t role in 



the formation of H2 (e.g.. 



Williams 



2OO3I) and H2O (e.g.. 



HoUenbach et al. 



2OO9I). Simple 
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theoretical arguments suggest that, even over the range covered by the Milky 



MCs, coohng and chemistry rates and timescales are affected by metalhcity (IBanerji et al 



Way and the 



20091 ). Not only is it crucial to understand the ISM chemistry at low metallicity, it is also 



important to investigate the chemical evolution during the star formation process. 

In cold molecular clouds, layers of ice form on the surface of dust grains. In particular, 
when densities increase dramatically during the cold collapse phase, most molecules accrete 
onto dust grain surfaces. In fact, as water, CO, and also CO2 are completely frozen out of 
the gas phase in these envi ronments, ices are subs tantial reservoirs of the available oxygen 



and carbon (see review by 



Bergin &: Tafalla 



20071 ). Even if CO seems to form exclusively 



in the gas-phase, there is lab oratory evidence th at water ice is predominantly formed by 



hydrogenation on icy grains ( lloppolo et al. 



20081 ). Therefore, understanding ice chemistry 



is crucial to understanding the gas-phase chemistry and the physical conditions within 
molecular clouds. Water is by far the most abundant ice (typically 10~^ — 10^^ with respect 
to II2), follo wed by CO? and C O, with a combined abundance of 10—30% with respect 



to water ice (IGibb et al. 



20041 ). Other ice species have also been identified (methanol. 



ammonia, ammoniu m, methane and formaldehyde), albeit in much smaller concentrations 



(Ivan Dishoeck 



20041 ). While water and CO can be observed from the ground, CO2 can only 



be observed from space and its ubiquitous presence in a variety of environments was one o f 



themajor discoveries of the Infrared Space Observatory (ISO) (see 



2003 



Ehrenfreund &: Eraser 



, and references therein). 



It is still unknown whether metallicity affects ice chemistry. The most obvious effect 
of a lower metal content is a reduced availability of carbon and oxygen, the cornerstones of 
interstellar chemistry. Indeed the CO abunda nce in the diffuse ISM in the LMC is lower 



than that measured in Galactic environments (Eukui et al 



2OOII ). In metal-poor clouds, it 



could be expected that ice abundances are generally depressed (i.e. all ices could be less 
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abundant as a result of lower abundances of C and O), but it is possible that the abundance 
balance of ice species is also reshaped. The lower dustigas ratio should make cloud cores 
more transparent to ambient UV light, possibly exposing ices to more radiative processing, 
and making complex species more abundant than in the Galaxy. 

Interstellar ices have been detected in the envelopes of heavily embedded young stellar 
objects (YSOs) in the MCs. The first evidence of ices in the env elope of a massive YSO in 



any extra-galactic environment was a serendipitous discovery by 



van Loon et al 



fl2005ah . 



The s pectrum of IRAS 5328— 6827 in the LMC, obtained with the Infr ared Spectrorn eter 



(IRS 



Houck et al 



20041 ) onboard the Spitzer Space Telescope {Spitzer, 



Werner et al 



20041 ) and the Infrared Spectrometer And Array Camera at the Very Large Telescope 



(ISAAC/ VLT), shows c l ear ab sorpt ion signatures of water ice at 3.1 /im and CO2 ice at 



15.2 fim. 



Oliveira et al. 



(120061) and 



Shimonishi et al 



(I2OO8I ) identified more embedded 



YSOs with ic e signatures in the LM C. Three embedded YSOs have been recently identified 



in the SMC 



van Loon et al. 



field (IWeltv et al 



20081 ). There is a suggestion that the stroriger U V radiation 



20061 ) and/or higher dust temperature (lAguirre et al. 



20031) could be 



responsible for a higher C02-to-water colum n density ratio o 



3served towards a few YSOs in 



the LMC when compared to Galactic YSOs (Ivan Loon et al. 



2005a; 



Shimonishi et al. 



20081). 



These studies are however hampered by small sample size and limited spectral resolution 
and thus not yet conclusive. 

CO2 ice is found to be abundant and ubiquitous in_a variety of molecular environments 



in the Galaxy, from q uiescent clouds (Whittet et al. 



low-luminosity YSOs (IGerakines et al. 



I999I : 



20071) to the env elopes of high- and 



Pontoppidan et al. 



20081 ) . This molecule is 



generally believed to form exclusively via surface reactions, even thoug h its formation 
mech anism is not fully understood (for a recent discussion of this issue see ipontoppidan et al. 
20081 ). The most prominent CO2 solid state features are the 4.3-/im stretching and 15.2-/im 
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bending modes (e.g.. lEhrenfreund et al.lll997l ). Comparisons of observed 15.2-yum profiles 
witfi laboratory spectra suggest that CO2 i ces are dominated by components that are eithe r 



water- or CO-rich, rather than pure CO2 (iGerakines et al. 



1999 



Whittet et al. 



Pontoppidan et al. 



2008 



2OO9I ). For some YSOs , narrow substructure in the observed pro files provides 



evidence of thermal processing (e.g. 



Gerakines et al 



1999 



White et al. 



2OO9I ). Therefore 



the CO2 ice profile is a sensitive probe of its molecular environment. 

Another ice band in the IRS range (5— 7/im) is less well understood, with contributions 
not only from water ice at 6 fim but also amm onia (NH.s), methano l (CH3OH), formic 



acid (HCOOH) and formaldehyde (CH2O) ices (IBoogert et al 



20081 ). This superposition 



of multiple ice components is possibly responsible for the discrepancy between the column 
densities determined from the 3-yum (O— H stretching mode) and the 6-/im (O— H bending 



mode) water ice features, with the 6-/xm column density systematically higher (iGibb et al. 



2004 



Boogert et al. 



(NHJ) flBoogert et al. 



2008h. There is an ice feature at 6.8 fim attributed to ammonium 



20081 ) . This band also falls in a wavelength range strongly affected by 



emission features generally attributed to polycyclic aromatic hydrocarbon (PAH) molecules. 
The libration mode of water ice is a very broad feature that peaks at 13 fim and extends 
between 10— 30/im; in practice this feature is very difficult to isolate both because of its 
width and the fact that it is hard to disentangle from the more prominent silicate dust 
features. 

In a series of papers we will investigate ice chemistry in the envelopes of massive 
embedded YSOs in the low metallicity e nvironment of the MCs. The Spitzer Legacy survey 



of the LMC (SAGE, 



Meixner et al. 



20061 ) ide ntified large nurnbers of YSO candidates for the 



20081 ). The follow-up Legacy 



first time in any extra-galactic environment (IWhitney et al.l 
program with the Spitzer-IKS includes a variety of point and extended sources among which 



YSOs feature prominently (SAGE-Spec, 



Kemper et al 



20091 ). We selected the LMC sample 
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described here from both the SAGE-Spec program and the Spitzer archive, as those objects 
that have spectral energy distributions (SEDs) consistent with early embedded YSOs and 
exhibiting in their spectrum the absorption feature associated with CO2 ice, the most 
prominent ice band in the Spitzer-IRS spectral range. 

In this paper we present the first detailed analysis of CO2 ice profiles for massive 
YSOs in the LMC The paper is organized as follows. First we introduce the sample and 
discuss the issues of continuum determination and constraining the PAH spectrum. We 
then describe the observed CO2 profiles in terms of laboratory components and measure 
column densities. We also investigate features in the 5— 7/im range. Finally we interpret 
our results in terms of possible effects of low metallicity. 



2. Observations and data reduction 

The Spitzer-IKS spectra of our target sample were obtained mostly as part of the 
Spitzer SAGE-Spec Legacy program (PID: 4015 9). For three targets , IRS spectra were 



Kemper et al. 



(I2OO9I ) describe in detail 



recovered from the archive (PIDs: 3591, 3505). 
the original target selection, observing strategy, and data reduction for the SAGE-Spec 
program; we provide here only a brief summary. 

The IRS observations of the point sources in SAGE-Spec were performed in staring 
mode, using the Short-Low and Long-Low modules (SL and LL respectively). The 8- and 
24-/im fluxes were used to set exposure times, aiming at signal-to-noise ratios of ~ 60 in SL 
and ~30 in LL. Archival spectra mentioned previously were also obtained in staring mode 
with SL and LL. 

The spectra were reduced following standard reduction techniques. Flatfielded images 
from the standard Spitzer reduction pipehne were background subtracted and cleaned of 
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"rogue" pixels and artifacts. Spectra were extracted individually for each data collection 
event (DCE) and co-added to produce one spectrum per nod position. The nods were 
combined to produce a single spectrum per order, rejecting "spikes" that appear in only 
one of the nod positions. Finally, the spectra of all the segments were combined including 
the two bonus orders, that are useful in correcting for discontinuities between the orders. 
Archival spectra of LMC objects were re-reduced by the SAGE-Spec team to the same 
uniform standard as the SAGE-Spec sources. 

When comparing spectra obtained with different instruments, and of objects at 
different distances, it is important to consider potential differences in the physical scales 
that are sampled. The Spitzer IRS shts are 3.6" and 10.5" wide, respectively for the SL 
and LL modes. In the cross-dispersion direction the typical 2-pixel scales sampled are 3.6" 
and 10.2", respectively. At the distance of the LMC {d ~50kpc), these approximately 
square regions on the sky correspond to 0.9 x 0.9 pc^ and 2.5 x 2.5 pc^ at 6 — 7 and 15 /xm, 
respectively. 

By comparison, the ISO SWS apertures are 14" x 20" around 3 and 6 — 7//m, and 
14" X 27" around 15 //m. Observations obtained with these modules of massive-star 
forming regions in the Galaxy at typical distances of ~ 1.5 kpc sample physical scales of 
0.10 X 0.15 pc^ and 0.10 x 0.20 pc^ at 6—7 and 15 /xm, respectively. Hence these SWS 
observations would sample scales roughly ten times smaller in typical Galactic regions than 
the Spitzer observations in the LMC. Observations obtained with Spitzer of low- mass YSOs 
in nearby complexes at typical distances d ~150— 350pc sample even finer scales, by a 
factor ~ 200 compared to observations in the LMC. 
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3. The massive YSO sample 



This study focuses on the least evolved, most embedded massive YSOs in the LMC 
Embedded YSOs can be selected based on their mid-IR colors. However, the colors and 
even broad-band SEDs of YSOs are in some instances simil a r to t hose of evolved objects 
like planetary nebulae and post-AGB stars. 



Whitney et al. 



(120081 ) discuss in detail the 



YSO selection criteria using Spitzer data and contamination by evolved stars. One of the 
main goals of the SAGE-Spec program is to use spectra to clarify the boundaries between 
different classes of objects in the color-color and color-magnitude space. For this study 
we have selected only objects whose nature is spectroscopically unambiguous, which also 
means that these objects are the most embedded and luminous. 

We initially identified YSO candidates based on their very red continuum over the 
5— 38/im range. Then we narrowed the sample to include only those objects that showed 
an indication of the CO2 ice fe ature at i5.2Atm. CO9 in the gas phase has been detected 



towards many AGB stars (e.g.. 



Justtanont et al 



19981 ): CO2 ice however is only observed 



in YSO envelopes and quiescent molec ular clouds (Section 4 .3) — water ice is also observed 



in the envelopes of evolved stars (e.g. 



Dijkstra et al 



2OO3I ). Therefore the shape of the 



SED and the CO2 ice feature unambiguously identify these objects as embedded YSOs. We 
did not make use of the 5— 7/im ice complex for sample selection, because this wavelength 
range is dominated by PAH emission in our YSO sample (see below). 



which were c 



assified as "high- 



20081). TablesEand 



The final sample then includes 15 objects, 10 o 
probability YSOs" based only on their IRAC colors flWhitney et al. 
[2] provide the positions and observed fluxes for the YSOs. Even if some objects are 
known IRAS or MSX sources, many objects had not been classifled or studied prior 
to the Spitzer surveys. Van Loon et al. (2005a) analyzed the IRS and ground-based 
L-band spectra of IRAS 05328— 6827, which show ice signatures; we re-analyzed here its 
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IRS spect rum 



(3—5 /iin, 



Shimonishi et a. 



Murakami et al. 



(120081 ) obtained low-resolution AKARI near-IR spectra 



20071 ) of three of the sources, which also showed ice features. 



Van Loon et al. (2005 b) tentatively identifie d IRAS 05246-7137 as an AGB or post-AGB 



object. More recently, 



Whitney et al. 



(120081 ) classified this source as a YSO based on its 



IRAC colors. The object's IRS spectrum (showing strong silicate absorption, weak PAH 
emission, and a tentative ice detection) and IR luminosity support the YSO classificati on. 



MSXLMC464 was tentatively identified as a possible Hll region (IKastner et al. 



20081) 



based o n 2MASS and MSX colors. Van Loon et al. (2009) describe the Spitzer MIPS-SED 



spectra (IRieke et al. 



200J) of IRAS 05328-6827, MSX LMC 1786, and IRAS 05246-7137; 



they all show a cold dust continuum between 53—92 /im consistent with a YSO status. 



4. Analysis of the IRS spectra and ice absorption features 

Figured] shows the spectra of the 15 objects analyzed in this work, covering the whole 
IRS range. The insets show the spectral regions associated with the ice bands, CO2 at 
15.2 /im and the 5— 7/im complex. Objects are displayed in decreasing fiux intensity at 
5.8 /im. Many objects show broad absorption at 10 /im (and occasionally 18 /im) attributed 
to silicate dust grains; one object (IRAS 05240— 6809) shows sihcate emission. PAH 
emission is also present in many objects in our sample: the prominent 6.2- and 11.3-/im 
PAH features are easily identified in at least two thirds of the spectra. Figured] also shows 
that there are a variety of strength and profile shapes associated with the CO2 ice and that 
few objects exhibit clear ice absorption at 6 /im. In the next sections we will analyze these 
features in more detail, starting with the issue of continuum determination. 



- 11 - 



4.1. YSO properties: SED fitting 

We have fitted the observed SE Ds of the LMC YSOs using the on-line fitting too 
developed by 



Robitaille et al 



13 



(120071 ) ■ This tool relies on a database of 20,000 pre-calculated 
models that are compare d to the observed SEDs, and the best fit is selected using a 



min imization te c hniqu e . 



and 



Chen et al. 



Robitaille et al. 



20071 ) explain these tools in more detail. 



(120061 . 

(j2009l ) give examples of their use. A typical YSO model includes a central 
stellar photosphere, a flared accretion disk, and a rotationally flattened infalling envelope 
with bipolar cavities, characterized by a total of 14 parameters. The balance of the different 
model components determines the object's evolutionary status. 

To constrain the observed SEDs we made use of the IRS spectrum discussed here as 
well as all the photometr ic data available for e ach object (Tables[T] & [2]). The SEDs of 10 



objects are described by 



Whitney et al. 



(120081 ): for the remaining objects the results of 
the SED fits will be discussed in detail by Chen et al. (in preparation). The objects in 
the sample have estimated YSO luminosities in the range ~ 5 — 50 x 10'^ Lq and masses 
between ~ 10 and 25 M0. Taking into account their observed luminosity, the objects in our 



sample are the 



Gerakines et al. 



VIC c ount erparts of the luna inous Galactic YSOs such as those studied by 



(Il999h and 



Gibb et al. 



(120041 ) with ISO. In terms of evolutionary stages the 



YSOs are all in Stage I, i.e they are still deeply embedded in their envelopes, as expected 



from the presence of ice absorp 
the Class I as defined by 



tion 



eatures. The Stage I classification is equivalent to 
Ladal (jl987l ). but instead of relying on the spectral index of the 
SED, it is based on modelled SED quantities, in this case the accretion rate of the envelope 
normalized to the stellar mass, Menv/M^ > 10^6 yr~i (lR,obitaille et al.ll2006[ l. 



^SED fitter available at |http:/ /caravan. astro. wisc.edu/protostars/index.php 
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4.2. Dust continuum, PAH emission and silicate features 

The mid-IR spectra of embedded YSOs exhibit a wealth of solid-state and gas features 
super-imposed on a dust continuum: silicate absorption/emission features at 10 and 18 /im, 
PAH emission features from 5 to 19 /im, molecular hydrogen and forbidden line emission 
across 5 — 35/im and finally the ice signatures we are interested in here. In our sample, 11 
out of 15 objects show prominent PAH emission (look for the conspicuous 11.3-/im feature 
in Figured]). As we want to measure ice optical depth, we need both to determine the 
continuum level and account for any emission or absorption features that fall in the desired 
spectral ranges. We fitted continua between 5 and 20 /xm, focusing on obtaining good 
continuum determinations around 5 — 7 and 15 /xm. 

Determining a global IR continuum for YSOs is not trivial. This is because the lines 
and bands listed above merge together, meaning that few true continuum points exist across 
the spectrum. As a consequence continuum fits rarely offer unique solutions. The CO2 ice 
band at 15.2 yum is relatively narrow and falls in the region between the 11.3— 12.7 /im and 
17-/im PAH complexes and between the silicate absorption bands. Therefore, it is possible 
and practical to use a local continuum determination, for instance a low-order polynomial 
fitted to the local pseudo-continuum in the neighborhood of the ice feature. 

This is not the case for the 5— 7/im ice band. First, there is prominent PAH emission 
in this region, namely the 6.3-/im feature and the 7.7-M m complex, a s well as unresolved 



Spoon et al 



fl2002[ ) show that 



H2 emission at 5.51 (S7), 6.11 (S6), and 6.91 /im (S5). 
relatively weak water ice signatures at 6 /im are difficult to separate from the PAH emission 
contribution. Secondly, this range falls at the edge of the IRS spectrum, where the spectra 
are noisier. This means that it is very difficult if not impossible to find bona-fide continuum 
points to anchor a local pseudo-continuum determination. 



PAHFIT ( ISmith et al.l 120071 ) provides a convenient way of constraining the dust 



- 13 - 



continuum and extinction as well as accounting for PAH and gas emission. The dust 
continuum is modelled with up to eight thermal dust components represented by modified 
blackbodies of a range of temperatures; dust extinction is produced by adopting extinction 
profiles similar to those of Milky Way dust. PAH emission features are represented by 
individual and blended Drude profiles — depending on whether PAH complexes are resolved 
into their multiple components — with fixed widths and central wavelengths. Molecular 
and atomic emission lines (e.g., H2 pure rotational lines, forbidden line emission from [S 
III], [Ne 11], and [Ne ill] etc) are modelled as Gaussian profiles; their widths and central 
wavelengths are allowed to change to account for the fact that these lines are not resolved 
in the IRS low-resolution spectra. To help PAHFIT constrain the continuum we added the 
IRAC 3.6- and 4.5-/im photometric points to the IRS range. The very broad water libration 
mode at 13 /im is treated as a component of the continuum in this approach. Once all the 
PAHFIT components are accounted for, we are left only with the ice contributions. 

Figure[2] shows fits to the 11 YSOs in the sample with PAH emission, using PAHFIT. 
Total fits are shown, as well as color-coded individual contributions of dust continuum, 
PAH emission and molecular and atomic line emission. For each object we show their 
spectrum and the continuum-subtracted spectrum to highlight the PAH contribution. We 
found that PAHFIT fits the spectra well from 5— 17yum. PAH emission is unequivocal in 
Figure[2l In terms of general appearance, the continuum-subtracted spectra of objects like 
MSXLMC464, SAGEIC J050354.56-671848.9, and IRAS 05452-6924 are reminiscent of 
the spectral comp onents betwee n 5 and 14 /im attributed to neutral PAHs, PAH clusters. 



and PAH cations (ITielens 



20081). 



At longer wavelengths, we find that PAHFIT has some trouble fitting the red side 
of the 17-yum PAH complex that also includes H2 emission at 17.0 fim (examples are 
IRAS 04514-6931 and IRAS F04532-6709 in FigureE]). In particular at about 18 /im there 
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seems to be more emission than can be produced by PAHFIT: there are no known PAH 
emission features between 17.9 and 18.9 fim. It is also not due to sp urious "features" relat ed 



Buchanan et al. 



2001 



to diffuse gas in the LMC (i.e. [S ill] at 18.7yum, as mentioned by 
The 17-/im complex is not as well studied as other, shorter-wavelength PAH features. Even 
though identified with ISO-S WS, it was on ly Spitzer-IRS that allowed a more detailed 



analysis of this complex (see 



Tielens 



2008 



for a review). It is described broad 



plateau that can extend from 16 to 19 ^m, but observations show a varied substructure. 



Van Kerckhoven et al 



(120001 ) and 



Peeters et al. 



(120041 ) have used large molecule databases 



to show that a huge range of merged profile structures can be achieved by combining 
individual PAH emission contributions. Thus it is not surprising that the fits are not perfect 
in this region. None of this affects the continuum underlying the ice bands. 

We note that a detailed analysis of PAH emission is beyond the scope of this paper, but 
we verified that PAHFIT is an adequate tool to use in the context of a YSO study. While 
the spectral contributions underlying the CO2 ice feature are well constrained, that is not 
the case for the 5— 7/im ice complex. Even using PAHFIT to account for emission features, 
it is very difficult to isolate this ice band: the only objects for which this was possible have 
no detectable PAH emission (see below). 

The four remaining objects do not show obvious PAH emission in their spectrum; we 
thus adopt the traditional approach of fitting a low-order polynomial to the continuum. 
Continuum points are selected to avoid the ice bands at 5 — 7 and 15.2 fim and the silicate 
feature at 10 fim. The 18-/im silicate feature and the water libration mode are very 
broad, and we do not try to extract them from the spectra; instead we treat them as 
contributing to the local "pseudo-continuum" . Figure[3] shows the continuum fits for these 
objects. The two objects at the top (IRAS 05246-7137 and IRAS 05328-6827) show strong 
silicate in absorption and weak CO2 ice absorption. SAGEIC J052546. 49— 661411.3 and 
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IRAS 05240—6809 show silicate emission and self- absorption respectively, even though the 
nature of the silicates is difficult to constrain without a detailed study. In general the 
silicate optical depth is not well constrained as the shape of the continuum underlying 
it is somewhat uncertain. Both these objects show ice absorption at 15.2 and 6 fim. 
IRAS 05240— 6809 also shows H2 emission at 5.51, 6.91, and 9.66 /xm. These two apparently 
PAH-free objects are the only YSOs for which we were able to successfully isolate the 
5 — 7 fim ice band. 

For the 11 PAH-rich YSOs we use the PAHFIT output to extract the optical depth 
of the ice features, for the remaining objects we use the spectrum after subtracting a 
polynomial continuum. We compared the optical depth spectra obtained when using the 
PAHFIT and the polynomial continua, for these four objects. This allowed us to constrain 
how the choice of continuum affects the measured column densities for the ice features 
(Section 4.3.3) 



4.3. CO2 ice at 15.2 /xm 



It has been shown in laboratory experiments that the CO2 ice profile at 15.2 /im 
is a sensitive probe of the environment in which the ice resi des, both in terms of 
physical conditions like temperature (it traces strong heating 



Ehrenfreund et al 



1998 



Gerakines et al. 



1999 



White et al 



20091 ) and also in term s of molecular environm ent (it 



Ehrenfreund et al. 



19971 ) as is 



can indirectly pinpoint other molecules in the ice matrix, 
discussed in the next section. On one hand this makes CO2 an interesting ice species to 
study, but on the other it implies that laboratory studies of this solid are very complex as 
the whole molecular system needs to be considered, not just CO2 ice in isolation. Realistic 
astrophysical environments are also impossible to replicate. Thus we cannot model the CO2 



profiles directly using experiments. The next best thing is a phenomenological approach 
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that attempts to model the observed profiles as a combination of "components" that have 
been well studied in the laboratory. 

CO2 does not form in any siz able abundance in the gas phase in the diffuse interstellar 



medium or molecular clouds (e.g., 



Herbst fc Leung 



19861 ). but it is not yet fully understood 



which of a number of possible grain-surface reactions is responsible for the abundance 
and ubiquity of CO';i ice s (for a recent discussion on CO2 formation mechanisms see 



Pontoppidan et al 



20081 ). It was initially thought that this molecn 



stron g UV photolysis of mixed ices containing H2O and CO (e.g., 



e would form durin g 



D'Hendecourt et al. 



19861 ). conditions easily found in the envelopes of lumin ous embedded YSOs. Ho wever, 



2007 



20091 ) indicate 



observations of abundant CO2 also in the quiescent ISM (IWhittet et al.l 
that additional radiative processing, besides that provided by cosmic-ray induced photons 
and ambient radiation field photons, is not required for CO2 formation. 



H2O ice is the most abundant ice species, followed by CO and CO2 (iGibb et al 



20041 ). These distinct ice species are often observed in the same sightline and are present 
on the same grain mantles possibly in a layered, "onion"-like structure. In particular, 
two independent components have been identified on icy mantles: H20-dominated and 
CO-dominated components, also referred to as p olar and apolar components, respectively 



(iTielens et al. 



1991 



Ehrenfreund &: Frasei 



2OO3I ) . A possible scenario is that H2O ice (with 



traces of CO, NII3, and CH4) forms easily and widely in molecular clouds, as it requires 
relatively low gas density and only moderately low temperature. At higher densities grains 
accrete mainly apolar molecules like CO (and IR inactive molecules like O2 and N2, see 
below): this rapid freezeout creates a mantle of almost pure CO on top of a water-rich layer 



(ITielens et al. 



1991 



Ehrenfreund et al. 



19971 ). 



Laboratory experiments show that 15.2-yum CO2 pr ofiles in ice matrices dom inated 



by either H2O or CO have very distinct properties (e.g. 



Ehrenfreund et al 



1997, see also 
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FigureHj). The comparison of laboratory and astronomical CO2 ice profiles shows that 



2008 



White et al. 



2009 



Whittet et al. 



by a combination 


polar (H^O-domina 


ted) and apolar 


Gerakines et al. 


1999 


Whittet et al. 


2007; 


Pontoppidan et al. 



20091 ). The polar compon ent usually dom i nates 



accounting for a large fraction of the total C O^^ column density ( iGerakines et al. 



1999 



Pontoppidan et al. 



2008 



Whittet et al. 



20091 ). Thus, even though doubts remain about 
its formation mechanism, these observations indicate that viable CO? format ion routes 



must be present both in H2O- and CO-dominated ice mantles (jPontoppidan 



2006 



Pontoppidan et al. 



20081). 



For some Galactic YSOs, the CO2 ice profiles show additional su 



3structure. The fact 



that this structure is absent in the CO2 profiles in quiescent regions (IWhittet et al. 



2007, 



20091 ) suggests th at thermal processing cou 



CO2 component. 



Ehrenfreund et al 



(11998 



d play a role in the production of this additional 



19991 ) showed that when heated the spectra 
of laboratory mixtures of H2O, CH3OH (methanol), and CO2 i n equal parts deve lop a 



complex structure due to annealing and crystallization (see also iWhite et al.ll2009l ). The 



samples become very inhomogeneous and at temperatures of typically 100 K CH3OH and 
CO2 segregate and inclusions of pure CO2 are produced. As a r esult the profiles develop the 



double-peaked structure that is characteristic of pure-C02 ice ( lEhrenfreund et al 



19971). A 



red shoulder at 15.3 fim also appears, very likely due to the interaction of CO2 and CH3OH 



molecules. Several studies have found that the annealed ice component combined with t 



he 



dominant water-rich component fits the obs erved CO2 ice profiles well (jGerakines et al. 



1999 



White et al. 



2009 



Zasowski et al. 



20091 ). Recently 



Pontoppidan et al. 



(120081 ) suggested 



that for the case of low-luminosity YSOs, where the relatively high temperature needed 
for the segregation process may not be reached, the pure-C02 component arises from the 
distillation of the C0:C02 component through evaporation of CO at lower temperatures 
(<50K). 
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Bringing all this i nformation together, we have the following schematic scenario 



(iGerakines et al. 



1999 



Pontoppidan et al 



20081 ). CO2 is formed plentifully on icy grain 



surfaces in quiescent molecular clouds, mainly in an H20-rich medium but with a 
contribution of a CO-rich mantle in agreement with observations. This complex ice mantle 
is heated up by the contracting YSO, and what happens next depends on the temperature 
profile of the dusty envelope. In the environments of massive luminous YSOs, high enough 
temperatures can be reached (~100K in laboratory experiments, 50— 80K in collapsing 
dusty envelopes) so that annealing and segregation gi ve rise to co-existing H-^O-rich and 



C02-rich (~ equal parts H2O and CO2) components ( IGerakines et al.lll999l ). This heating 



process would sublimate the CO-rich componento. In the envelopes of lower luminosity 
YSOs, only moderate he ating takes place (5, 50K) ; the CO gradually desorbs leaving behind 



a pure-C02 component (jPontoppidan et al. 



20081 1. In the latter case, H20-rich, CO-rich 



and pure-C02 components co-exist in icy mantles. How these two scenarios are played 
out in a YSO environment would depend on the luminosity and evolutionary stage of the 
object, and the temperature profile across the dusty envelope. 



Recent st udies have favored com bining the water-rich component exc 



annealed ices (IZasowski et al. 



usively with either 



20091 ) or CO- and/or C02-dominated ices (jPontoppidan et al. 
20081 ) to model the observed profiles. Both approaches fit the data well, and it is not clear 
which of these is more appropriate. We avoid making a priori decisions on the nature of the 
profiles observed towards the LMC sources. We model the CO2 ice profiles in our sample 
in both ways: i) by combining H20-rich, CO-rich, and pure-C02 components; and ii) by 
combini ng a H^O-rich compone nt w ith an annea l ed com ponent — the same approach is 



used by 



Gerakines et al. 



(11999h and 



White et al. 



(120091 ). The next sections describe the 



laboratory ice profiles used in the analysis and the resulting profile fits. 



^Sublimation temperatures are of the order of 20—30, 50—90, and 100— 120 K, respec- 
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4-.3.1. Laboratory ices 



The laborator y ice p rofiles used in this study come from public databases from 

Leiden database: polar, apolar, and pure CO2 ices) and 



Ehrenfreund et al 



White et al. 



(119971 



(I2OO9I . annealed ices), which Tabled lists. In the Leiden database. 



particle shape correction s are applied to the labor atory profiles using di: 



(Ehrenfreund et al. 



19971 ) ; following other authors (iGerakines et al. 



1999 



erent grain models 



Pontoppidan et al. 



20081 ) we use a Continuous Distribution of Ellipsoids (CDE, ea ch grain shape equally 



probable) that seems to work best for CO2 profiles. According to 



Gerakines et al 



(Il999h no 



such particle shape corrections are appropriate for the highly inhomogeneous annealed ices. 

These databases are incredibly rich, and many parameters can be tuned (laboratory 
temperatures, fractions of ices in mixtures, etc). Instead of using the full da tabases and just 
testin g which combination of lab profiles fits each spectrum best (as done in 



Gerakines et al. 



19991 ). we simplify the problem (jPontoppidan et al.ll2008l ). and reduce whenever possible 



the number of free variables — i.e. temperatures are kept fixed except in the case of the 
processed ices where thermal effects are being probed, and relative concentrations are 
only changed when it affects the morphology of the profiles in a unique way. This should 
make possible trends easier to detect. The lower resolution of the Spitzer LL mode (when 
compared to previously published analysis) also justifies this simplification, as profile details 
are lost (see below). F or the same reason we on ly adopt the three dominant components 



used in the analysis of 



Pontoppidan et al. 



(120081) (HaO-rich, CO-rich, and pure CO2). In 



the rest of this section we describe the profiles used. 

For polar (water-dominated) ices, two mixtures are available in the databases: 



tively for CO, CO9, and H^O depending on 



(iGerakines et al. 



199S 



van Broekhuizen et al 



:he ice mixtures in whi ch the species reside 



2006 



Eraser et al. 



200 ll ) 
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H20:C02 = 100:14 at 10 K and H20:C02:CO = 100:20:3 at 20 K. As it can be seen from 
FigureH] (top left) these two mixtures produce very similar spectra, with only a minor 
shift in the peak position; still we consider them both for completeness sake. A 4*^-order 
polynomial is used to remove the contribution of the H2O libration mode at 13 /im from the 
spectra of all water-rich mixtures used. 

For apolar mixtures of CO and CO2 several types of mixtures are available, either 
just C0:C02 or including IR inactive molecules, CO:02:C02 and CO:N2:C02. IR inactive 
molecules like O2 and N2 are formed in the gas phase and accrete onto the cold grains at 
higher densities. They have no innate IR transitions, but in an ice matrix they interact w ith 



19971 ). 



other molecules, subtly changing the shape of strong profiles ( lEhrenfreund et al. 
For the same C02:C0 concentration ratio, the presence of O2 in the ice matrix makes the 
CO2 ice profile slightly wider (FigureH] top right, top two pairs of spectra). However, the 
exact amount of O2 in the mixture is not well constrained (profiles from spectra where the 
amount of O2 is 10, 20 or 50% of CO are indistinguishable). Therefore we keep the O2 
concentration at 50% of that of CO. The presence of does not change th e profiles in a 



Ehrenfreund et al. 



19971 ) so we do not 



noticeable way (N2 is a "silent" matrix component, 
include those laboratory profiles. 

With increasing CO2 fraction in the ice mixtures, the CO2 profiles become progressively 
more asymmetric, with a more pronounced blue wing and also blue-shifting the peak 
position (FigureH] top right, bottom two pairs of spectra). These changes in the profile shape 
are monotonic with 0^:00 ratio; thus we can interpolate between available laboratory 



spectra and ratios (see 



Pontoppidan et al. 



20081 ) . Available spectra have C02:C0 ratios as 



listed in TableO The temperature of the mixture (generally 10, 20 and 30 K) also has no 
unique effect on the shape of the profile for any of the mixtures mentioned so far; we adopt 
the profiles for 10 K. 
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Laboratory spectra of pure-C02 ices are available at three temperatures (10, 50, and 
80 K); with increasing temperature the profile becomes narrower, and the component at 
15.25 yum becomes more pronounced (FigureHI b ottom left). We us e the laboratory spectra 



White et al. 



(120091 ) as representative of 



for the H20:CH30H:C02~ 1:1:1 mixture from 
the annealed ice component. At lower temperature those profiles are rather broad and 
structureless; as the samples are gradually heated the characteristic peak splitting becomes 
pronounced (>100K), and the profile becomes considerably narrower (FigureH] bottom 
right). Profiles are available at temperatures from 5 to 135 K in 5-K steps (Table[3]). We 
should point out that laboratory temperatures quoted throughout this paper correspond 



to lower temperatures in conditions typical of dense molecular clouds (lEhrenfreund et al. 



1998 



Boogert et al. 



2000). 



A crucial thing to reahze is that previous studies made use of spectr a with higher 



spect ral resolving power. Observations per: 



1999) and t 



le Sy itzer IRS SH mode (e.g. 



White et al. 



brmed with the /5'0-SWS (e.g., 



Pontoppidan et al 



2008 



Gerakines et al. 



Zasowski et al 



2009 



20091 ) have 1000 and 600, respectively; for the IRS LL mode at 15.2 /xm 
90. We take this into account by convolving the laboratory ice profiles described 
previously with a Gaussian profile at each wavelength element with a width appropriate to 
simulate the lower resolution of the LL mode; FigureH] shows the effect of the resolution on 
the profiles. The convolved profiles for pure CO2 and annealed ices at higher temperature 
(lower two panels) show that the substructure in the profile that is normally associated 
with thermal processing is smoothed out at the LL resolution. The double-peaked structure 
is no longer clearly visible. Thus we are not able to diagnose processing in our spectra as 
readily and to the same level of detail as in previous studies. 

FigureH] also shows which particular components contribute to specific parts of 
the profile. For illustration purposes, we divide the profile into a blue and a red wing 
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(respectively short- and long-wards of 15.2 yum). The H20-rich and annealed mixtures at 
low temperature contribute to both wings of the profile and the red wing is much more 
extended than the blue wing (~ 1.1 versus 0.45 /xm). The CO-rich, pure-C02, and annealed 
mixtures at high temperature contribute almost exclusively to the blue wing of the profile 
— the pure CO2 profiles extend the furthest into the blue. 



4-3.2. Laboratory fits to the CO2 profiles 



The IRS spectra in our sample were compared to the laboratory profi les using a 



minimization technique (Levenberg-Marquardt least-squares mini mization 



More 



making use of the IDLo routines developed by 
discusses the resulting fits. 



Markwardt 



(12009 



19771 ) 



MPFITH). This section 



Each observed CO2 profile is modelled in two ways as described previously. Figure|5] 
shows the resulting fitted profiles. TableH] lists relevant fitted parameters as well as the 
total CO2 column density (see next section) and the percentages (in terms of column 
density) for each component. For three YSOs (MSXLMC464, IRAS 05246-7137, and 
SAGEIC J054059. 29— 704402.6) the optical depth spectra were not deemed of good enough 
quahty to fit meaningful models. Figure[6] shows their CO2 profiles, and TableH] lists the 
column densities. 

The spectra of the remaining 12 YSOs in the sample are well fitted by the 
combinations of the lab profiles. Some objects exhibit features at ~15.6/im (e.g., 
SAGEIC J052546.49-661444.3 and MSXLMC1786) and ~16.2/xm (IRAS 04514-6931) 
that none of the model components can reproduce (FigureHl). There is a [Nelll] emission 



^Interactive Data Language, 
"^http: / / purl.com/net / mpfit. 



- 23 - 



line at 15.56 yum and a weak PAH feature at 15.9 /xm. Therefore it is possible that these 
features are artifacts related to these two emission lines, either from the object itself or 
from its sky background. 

The fitted profiles show that we are not able to distinguish between the different H2O- 
and CO-rich mixtures used: similarly good fits are possible if the H20:C02 component is 
replaced by the H20:C02:CO mixture or the C0:C02 is replaced by CO:02:C02. There 
is also only a negligible difference in the fits when using the three available temperatures 
for pure-C02 ice. Accordingly, we make no such distinctions in TableHj referring only 
to H20-rich (polar), CO-rich ("apolarl") and pure-C02 ("apolar2") components. When 
the polar + annealed combination is used, laboratory temperatures > 100 K are required, 
consistent with thermally processed mixtures in which CO2 ice is segregated. 

For seven YSOs the two strategies provide equally good results in terms of 
minimization. For two YSOs the polar + apolar modeling fits their spectr a better while 
for th ree other objects the polar + annealed combination performs better. 



Gerakines et al. 



( 119991 ) found that the combination of polar + annealed components fitted the data of their 
Galactic sample better, and suggested that this indicated that ice processing has occurred. 
We believe that for our data we cannot make this particular discrimination. The laboratory 
components are obviously very different, but the loss of spectral detail at the LL resolution 
and the combination of multiple components limits what we can discriminate. 

In most cases, the two methods agree well on the amount of water-rich ice required. 
The majority of YSOs have a large fraction of CO2 locked into H20-rich ice: for eight of 
the twelve fitted objects the polar component accounts for ~45 — 75% in column density. 



This is consistent wit h what is observed for Galactic YSOs (e.g.. 



Gerakines et al. 



Pontoppidan et al. 



2008 



19991 ). For the other four objects either no water-rich ice is required or 



the two strategies do not agree on how much is necessary. The polar component provides 
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the dominant contribution to the red wing of the CO2 profile (FigureH]). When contrast in 
this region is low (as is the case for SAGEIC 045228.65-685451.3, IRAS 05452-6924, and 
IRAS 05328— 6827), it becomes more difficult to determine the contribution of the polar 
component. We test this effect by adding noise to one of our best profiles and re-fitting it: 
the profile looks narrower, and the amount of the polar CO2 ice is not well constrained. 
The spectrum of SAGEIC J052546. 49— 661444.3 is of good quality, but it is also dominated 
by the feature at ~ 15.52 fim. 

If we consider the polar + apolar fits, a contribution from pure-C02 ice is present for 
all but two fitted YSOs; it typically accounts for 20 — 50% of the t otal column density . 



Simila rly, the CO-rich component is present in 10 profiles. Unlike 



Pontoppidan et al. 



( l2008l ). we have difficulty in constraining the C02:C0 concentration ratio for the CO-rich 
component in the profiles individually. However, there seems to be a tendency in the 
sample for high C02:C0 ratios. If we consider the polar + annealed fits, a contribution from 
processed ices is needed in all cases. 



4-3.3. CO2 ice column density 



Table HI lists computed colum n densities adopting a band strength A = 1.1 x 10 ^''cm 



molecule ^ (IGerakines et al 



19951 ): CO2 column densities are in the range 0.45 — 16 x 10^^ cm~^. 
Uncertainties listed in the table are statistical only. Different choices in defining the under- 
lying continuum can change the measured column densities by ~ 15%. 

We can compare our new CO2 column density measurements for IRAS 05240— 6809, 
SAGEIC J052546.49-661411.3, and SAGEIC J051449.41-67221.5 with measurements 
performed with the lo w-resolution spectrogra ph (i?~20) on board AKARI for the 4.3-/im 



CO2 stretching mode (jShimonishi et al 



20081 ) . Their measurements suffered from very large 
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uncertainties (at least 50%, see TableH]); both sets of measurements are consistent within 
the quoted uncertainties, however our measurements bring the estimated CO2 column 
densities down considerably. 



For those three objects. 



Shimonishi et al. 



(I2OO8I ) also obtained water ice column densities 



— for IRAS 05240—6809 they do not tabulate a value for the column density, but we estimate 
it roughly from their Figure 1, ~ 40 x 10^'' cm~^ with an assumed error of 25%. Van Loon et 
al. (2005b) measured the column density of water ice in IRAS 05328— 6827; we re-estimate 
this column density using their spectrum and obtain 3.26 ± 0.27 x 10^'' cm^^. For these four 
YSOs (IRAS 05240-6809, SAGEIC J052546.49-661411.3, SAGEIC J051449.41-67221.5 
and IRAS 05328—6827) we combine our CO2 measurements with the H2O measurements. 



Shimonishi et al. 



(I2OO8I ) computed CO2 and water column densities for a further three 



objects. Thus there are now seven YSOs in the LMC with column densities measured for 
both ice species. There is a strong positive correlation between the CO2 and H2O column 
densities for Galactic YSOs; in the next section we compare the measurements for YSOs in 
the LMC and the Galaxy. 



4.4. The 5—7 fim ice band 



The 6- and 6.8-Atm ice bands dominate the 5— 7/im range in YSOs. Even though first 



identified by 



Puetter et al 



(I1979I ). their carriers are still not conclusively ide ntified. The 



mairi contributor to the 6-/im component is the bending mode of water (e.g.. 



Tielens et al 



20041 ). However, the column densities for the 3- and 6-/xm w ater bands were fo und to 



be discrepant, with the 6-fim band often much too deep (e.g. 



Keane et al. 



200 ll ). Other 



ice species like formic acid (HCOOH), formaldehy de (CH^O), and ammonia (NHQ may 



contribute to the de pth of the 6-/im feature (e.g.. 



Zasowski et al 



Knez et al 



2005 



Boogert et al. 



2008 



20091 ). Further, the 6-fim bending mode might be sensitive to the presence 
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Bouwman et al. 


2007 


), which 


ion at 6 /xm ( 


Knez et al. 


2005) 



The ice component at 6.8 fim is now beheved to 



3e due mostly 



a smaller contribution from methanol (CH3OH) (IBoogert et al.ll2008l ). 



to am monium (NH^ ), with 



In Section 4.1 we discussed how PAH and H2 emission have also important spectral 
signatures in the 5 —7 /im spectral range, in particular the strong 6.3-yum PAH feature. 



Spoon et al. 



(I2OO2I ) have shown that weak water ice features at 6 /im are effectively masked 
in the presence of a strong PAH emission spectrum (their Figure 1). As a large fraction 
of our sample shows significant PAH emission, this effect limits our ability t o identify ice 



absorption in this spectral range. A quick inspection of the LMC sample of 



Seale et al. 



( I2OO9I ) reveals that the majority of the objects that show CO2 ice also show prominent PAH 



emission. 



Boogert et al.l (120081 ) and IZasowski et al.l (120091 ) recently analyzed IRS high-resolution 
spectra of Galactic YSOs and found little or no sign of significant PAH emission. However 
these samples are mainly of low-luminosity objects. Thus PAH emission, if excited by the 



central ob 



e ct, wo uld not dominate. In the higher-luminosity ISO sample discussed by 



Gibb et al. 



(120041 ) some objects show PAH emission, even though not as prominent as in 



our LMC sample. 

In s pite of the fact that there is a global dearth of PAH emission in metal-poor 



galaxies ( IMadden et al. 



20061 ). we believe that the reason why PAH emission is so readily 
observed in YSOs in the LMC has to do with the issue of spatial resolution. As discussed in 
Section 2, observations YSOs in the LMC probe a larger region of the objects' environment 
compared to Galactic sources. Extended PAH emission can then easily appear superimposed 
onto the spectra of icy mantles. A Galactic example of such "source confusion" in star 
forming environments is MonR2: its ISO-SWS spectrum shows contributions from both 
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an embedded YSO with ice and silicate absor ption (MonR2IRSl) and an ultracompact 



Hll region (MonR2IRS2) fISpoon et al 



2004 . their Fig 4). Such problem is also well 



documented in the ana lysis of the spectra of ultra-luminous infrared galaxies (e.g., 



Spoon et al. 



2002 



20041 ). 



Figure[7] shows the spectra for the two objects in the sample for which we 
could identify the H2O ice absorption band at 6 yum: IRAS 05240 —6809 and 
SAGE1CJ052546. 49— 661411.3. None of these spectra show obvious PAH emission 
but IRAS 05240 —6809 shows unresolved H2 emission. These YSOs show two of the largest 
CO2 ice column densities (>9 x lO^^cm"^, TableS]), and thus should also have the largest 
water ice column densities (see next section). IRAS F04532— 6709 shows similarly large CO2 
column density, but we cannot detect the 6-/im water ice feature which is probably hidden 
by PAH emission. 



We also show in Figure[7] typical laboratory ice spectra (T = 40 K, 



Hudgins et al. 



19931 ) scaled to match published water column densities measured for the 3-/xm feature 
(TableH]) and taking into account the measurement uncertainties (shaded areas in the 
figure). We have not made any allowance for the variable contribution of other ice species 
to the water ice band as described above. It is immediately obvious that, in particular 
between 6 and 6.5 fim, the observed features are not deep enough to match the published 
water column densities. This suggests a deficiency in our continuu m determination in th is 



region. The figure also shows an example of a Galactic YSO from 



Zasowski et al. 



fl2009h 



for the complete 5 — 7 fim range, highlighting the contribution of the 6.8-/im ice band. The 



spectru m of SAGEIC 



to NH^ flBooeert et al. 



052546.49—661411.3 does indeed show this feature, likely attributed 



20081 ). 



- 28 - 



5. Discussion 

In this section we compare in more detail the properties of ices in the LMC YSOs to 
Galactic samples and interpret our results in terms of the star formation environment in 
the LMC. 



5.1. Observed properties of the CO2 ice profiles 

As we have seen in Section 4.4, a large fraction of CO2 ice is locked in a water-rich ice 
matrix, with a possible contribution from CO2 ice in a CO-rich layer. Figure[H] compares the 
fraction of CO2 ice in the water-rich component for the LMC sample and Galactic sample s 



from the literature: high-luminosity YSOs and background source 



Whittet et al. 



20091 ) and low- luminosity YSOs (jPontoppidan et al. 



Gerakines et al. 



1999 : 



20081 ). The column 



densities for this component are derived using the polar -|- apolar modelling approach 
(Section 4.3.2) for all samples. In the Galaxy, luminous YSOs and background sources 
display a larger fraction of CO2 ice in the water-rich component (~ 85%) when compared to 
low-luminosity sources (~70%). The typical fraction for the LMC objects is ~60%. This 
suggests that, even though both in the LMC and in the Galaxy the water- rich component 
is the major contributor to CO2 ice, its fraction might be smaller for the LMC sources. 
However, we note that the LMC sample is small, and systematic uncertainties could bring 
the samples into better agreement. Even though a similar profile decomposition was used 
to isolate the water-rich contribution, there are differences in the modelling approach (e.g., 
number and properties of the components considered). 



A contribution from pure-C02 ice is also needed to fit the majority of the profiles. 



^Background objects are bright sources that coincidentally sit behind the molecular cloud 
but are not associated with it in any way. They represent quiescent sightlines. 
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Whether this contribution arises from inclusions in annealed ices or a pure-C02 layer 
resulting from CO desorption is less obvious. If CO2 forms in either water- and/or CO-rich 
ice mixtures as has been widely suggested, then the presence of pure CO2 diagnoses thermal 
processing. Therefore, our analysis of the morphology (i.e. shape and composition) of 
the 15.2-/im CO2 ice profiles of a sample of embedded YSOs in the LMC reveals only 
marginal differences when compared to Galactic YSOs. In the next section we discuss the 
concentration of CO2 ice in relation to water ice. 



5.2. CO2 ice abundance in LMC 

Using compiled column- density measurements (Figure[9]), we recompute the Galactic 
A^(C02)/A^(Il20) ratio. Filled and open circles are for Galactic YSOs and background 
source s, respectively. Measur ements for LMC YSOs are shown as star s ymbols. The YS Os 



in the 



Gerakines et al 



fl200lh 



Pontoppidan et al 



objects. The 



Knez et al. 



19991) sam ple are luminous, while t 



lose in the 



fl2008h and 



(120051 ) and 



Zasowski et al 



Whittet et al. 



(12007 



Nummelin et al. 



(120091) sa mples are low-luminosity 



20091 ) samples are exclusively of 



background sources. 



Measurements of t 


le column density of CO 


mode ( 


Gerakines et al.l 


1999 


■ 


Knez et al. 


2005 


; I 



2009 



Whittet et al. 



Shimonishi et al. 



Pontoppidan et al 



2009! ) or the 4.3-/im stretching mode (iNummelin et al 



2008 



Zasowski et al 



2001 



2OO8I). Measu rements of II2O ice are for the 3-/xm stretching mode. 



1 

except for 1 Zasowski et al.1 (120091) . who used the 6-/im bending mode. They did not account 
for the contribution of other ice species to the 6-/im optical depth. Therefore, they have 



probably overestimated the column density o: 
samples, we need to scale them accordingly. 



' water ice. To comp are them with the other 



Boogert et al 



(I2OO8I ) compared the column 



densities of water ice at 3 and 6 fim, and they found that the 6-/im value is on average 
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higher than the 3-fim value by a factor of 1.74 (the r aedian of 49 sources; the standard 



Zasowski et aL 



(120091 ) appear in 



deviation is 0.62). The adjusted column densities of 
FigureO We have also converted their column densities for CO2 ice (computed using 
A = 1.5 X 10~^'^ cmmolecule"^ ) to use the same band strength as the other measurements 
discussed here (A = 1.1 x 10^^^ cm molecule"^). 

Figure[9] shows the strong correlation between A^(C02) and A^(H20). The median 
value for the number den sity ratio is iV(C02)/iV(H20) ~0.2 for luin i nous Galactic YS Os 



and background so urces (IGerakines et al 



1999 



Whittet et al. 



2007 



Knez et al 



2005 



Whittet et al. 



(INummelin et al 



spread. 



2009h. F o r low-luminosity Ga l actic YSOs, that ratio is N (CO^) /NjR^O) ~ 0.3 



2001 



Zasowski et al 



Boogert et al. 



2004 



Pontoppidan et al. 



20081 ) with a very large 



(120091) reported a ratio of 0.16; however if the correction factor 
derived in the previous paragraph is applied to the measured water column densities 
then that ratio becomes 0.28, fully consistent with the previous value. In fact, the 
right-hand panels (b and d) in Figure[9] show that, with this correction, the low- luminosity 
measurements are essentially indistinguishable for the different samples. Note that there 
is a larger scatter for the measurements of low-luminosity objects when compared to the 
high-luminosity and background sample. Therefore, Figured suggests a clear difference 
between the high -luminosity YSO and ba ckground samples, and the low-luminosity YSO 



samples (see also 



Pontoppidan et al. 



20081) 



Shimonishi et al. 



(I2OO8I ) estimated a ratio A^(C02)/A^(H20) ~ 0.45 ± 0.17 for five LMC 
YSOs. As mentioned previously, we revised some of Shimonishi's estimates downwards. 
For the seven LMC YSOs with both column- density measurements we estimate a median 
ratio A^(C02)/A^(H20) = 0.32 ± 0.15 (Icr). A single Shimonishi's measurement drives this 
large spread, it would otherwise be A^(C02)/iV(H20) = 0.32 ± 0.05 (la). Estimated total 
luminosities in our sample are 5 —50 x 10'^ Lq, comparable to the massive YSO sample in 
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the Galaxy. Figure[9] (left panels) shows that if the high-luminosity Galactic and LMC 
samples are compared the A^(C02)/iV(H20) ratio is higher in the LMC, by a factor 1.6. 
It is curious, however, that the ratio measured for the LMC objects is consistent with the 
low-luminosity Galactic sample. Note that column-density measurements do not depend on 
any modelling details but are prone to systematic uncertainties, namely in the continuum 
determination, that are generally difficult to quantify. 

One should investigate the homogeneity of the samples discussed here. The LMC 
sample is restricted to very luminous YSOs, and it should, in terms of luminosity at 
least, be comparable to the sample of luminous Galactic YSOs observed with ISO. Spitzer 
studies of star formation concentrated on nearby star-forming regions, molecular clouds 
and isolated globules that form predominantly intermediate- and low-luminosity objects. 
These low-luminosity Galactic samples include a variety of evolutionary stages from more 
embedded YSOs to Herbig AeBe, T Tauri and zero-age main-sequence stars. The different 
timescales involved in the formation of objects with different luminosity may leave an 
imprint on the observed ice chemistry. Furthermore, YSOs in the LMC should in general 
appear more evolved due to the smaller dust column density in LMC objects when compared 
to Galactic objects. 

Galactic and LMC measurements sample different spatial scales (Section 2), from a few 
parsecs in the LMC to milliparsecs in the Galaxy. If there is a gradient of ice abundances 
between the more immediate YSO cocoon and the wider molecular environment in which 



it is immersed (jPontoppidan 



20061 ). then the fact that ice properties are integrated over 



different spatial scales could affect global abundance ratios. 

Such issues should be considered, since models of molecular clouds and their chemical 
evolution during the star formation process sugg est that both ternporal and radial variations 



of molecular concentrations are important (e.g.. 



Lee et al. 



2004 



Rodgers &: Charnley 



2003 



- 32 - 



Garrod et al. 



2008 



van Weeren et al 



2009 



HoUenbach et al. 


2009 


)• 


Lee et al. 


(2004) 



singles out CO abundances (and the balance of its gas and solid phases) as particularly 
critical in determining the chemistry of star-forming cores. Even though they did not 
explicitly discuss CO2 formation, the evolution of the two molecules is obviously closely 
linked. The balance of water in the gas phase and in ice form is also controlled by time- 



and space-dependent effects ( iHoUenbach et al. 



20091 ). 



5.3. The star formation environment in the LMC 



The focus of this project is how the metal-poor star-formation environment in the 
LMC influences the chemistry of the envelopes that surround YSOs. Thus we now look at 
what is known about the environment in molecular clouds in the LMC. 

The structure of Giant Molecular Clouds (GMCs) in metal-poor galaxies including the 
LMC is different from that in the Galaxy. Observations sugges t that in me tal-poor galaxies 
there are large reservoirs of H2 that are not traced by CO (e.g.. Ilsraellll997l ). The magnitude 



of this H2-to-CO excess was found to depend strongly on both metallicity and strength of 



the radiation field. This has been interpreted as the result of se 



Israel 



1997 



ective photodissoci ation 



Bolatto et al. 



I999I): H2 



of CO in the outer parts of metal-poor GMCs (e.g., 
readily self-shields while CO is shielded from the photodissociating radiation mostly by 
dust, which is less abundant at low metallicities. In this case, CO emission would trace only 
the inner parts of metal- poor GMCs. If CO cores are smaller in the LMC, du e to the harshe r 



UV r adiation field (e.g.. 



Weltv et al. 



20061 ) and lower dust abundance (e.g.. 



Bernard et al. 



20081 ) then we could naively expect CO2 production to be reduced. A lower production 



rate of CO o could just result from the 



LMC (e.g. 



Bel et al. 



1986 



Fukui et al 



ower abundance o 



2001 



" gas-p hase CO observed in the 



Pineda et al 



20091 ). It could also result from 



a smaller density of dust grains that act as seeds for the surface chemistry. Both effects 
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would then tend to reduce the value of A^(C02)/iV(H20). However as described above this 
is not what we observed in the LMC. 

The formation of CO2 ice in a wide variety of environments (e.g., both quiescent 
sightlines and low- and high-luminosity YSOs) argues against the need for enhanced 
radiative processing. Still, laboratory experimen ts do show that UV irradiation enhances 



CO2 production (e.g.. 



D'Hendecourt et al 



19861 ). and models of diffusive surface chemist ry 



can p roduce high CO2 densities at slightly elevated temperatures (e.g.. 



Ruffle fc Herbst 



200lh . The LMC has a stronger ambient UV radiation fi eld when compared to the Galaxy 



(by a factor of 1—10, 



Bel et al. 


1986; 


Weltv et al. 


2006) 



20061 ) . As a result of the stronger 



radiation field, the average ISM du st temperature is also generally high er in the LM C 



(12<Td<35K, 



Bernard et al 



Reach et al 



19961 ) 



120081) than in the Galaxy (19 < < 25 K 
— dust temperatures in cloud cores are lower than in the diffuse ISM. Both these effects 
could increase CO2 production and increase A^(C02)/iV(H20), consistent with what is 
observed. 

We should also consider the possibility that the formation of water ice is somewhat 



inhibited in the L MC. The largest co 



is ~ 5 X 10^^ cm ^ ( Shimonishi et al. 



umn density of water ice measured so far in the LMC 



■ow the hig 



lest column densities 



2004 



Pontoppidan et al. 



20081 ). This is be 
meas ured towards Galactic YSOs, 10 —30 x 10^^ cm~^ (iGibb et al 
20081 ). Due to the small number of LMC objects investigated to date and possible selection 
effects this may not be a significant difference. Water ice seems to form rather easily 
in molecular clouds. However the Ay-threshold for water ice formation {Ay ~ 3 mag for 



typical Galactic molecu 



and dust temperature (iHoUenbach et al. 



ar clouds) is affected b y the strength of the incident radiation field 
2OO9I ). Still, it is difficult to see how water ice 



formation could be inhibited when CO2 ice formation is not. 



6. Summary and final remarks 



We analyzed Spitzer-IRS spectra of a sample of embedded young stellar objects in the 
LMC, obtained as part of the SAGE-Spcc Legacy program. We selected sources for this 
analysis based on the shape of their mid-IR continuum and the presence of the bending 
mode of CO2 ice at 15.2 //m; a total of 15 sources were included, with three CO2 ice 
detections considered marginal only (<3-(t). Based on modeUing of the objects' SEDs, we 
characterize the LMC YSOs as luminous and still in the early embedded stages of evolution. 

A large fraction of the spectra of the LMC YSOs show prominent PAH emission. This 
does not necessarily mean that the sources are more evolved, ultracompact Hll regions. 
The Spitzer-IRS spatial resolution combined with the distance to the LMC implies that our 
spectra sample a relatively large fraction of the YSOs' environment (a few square parsecs) 
when compared to Galactic YSOs, and thus PAH emission frequently appears superimposed 
onto the spectra of embedded YSOs. This makes the identification of the 5—7 /xm ice band 
difficult. We detect water ice at 6//m for two YSOs and the 6.8-//m feature (hkely due to 
ammonium) in a single object. 

It is usually believed that CO2 forms on icy mantles, in water-dominated (polar) 
and/or CO-dominated (apolar) ice matrices. In Galactic YSOs, the water-rich component 
is the most abundant. These components exhibit a distinct morphology allowing a 
multi-component analysis of the observed profiles. The double-peaked structure found in 
the CO2 ice profiles towards YSOs in the Galaxy is a signature of thermal processing, 
via which the CO2 segregates from the other ice constituents. There are two types of 
modelling usually found in the literature: a dominant water-rich component is combined 
with either one or more apolar (CO- or C02-dominated) components, or with an annealed 
ice component (water, CO2 and methanol in equal parts). Accordingly, we model the 
observed profiles of CO2 ice towards LMC YSOs in these two ways, using databases of 
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laboratory ice profiles available in the literature. We find that at the lower-resolution of the 
IRS LL mode, a considerable amount of detail is lost from the laboratory spectra. Thus our 
comparison of the two types of modelling is not conclusive. 

We find that the CO2 ice towards LMC YSOs is also dominated by the water-rich 
component, and they also require a contribution from pure CO2 ice. Thus in terms of 
general morphology and composition the CO2 ice towards LMC YSOs is similar to that 
observed towards Galactic YSOs. There are however two possible differences. First, there 
is a hint that the fraction of CO2 ice locked into a water-rich matrix is smaller in the LMC. 
Furthermore, the column density ratio A^(C02)/A^(H20), available for some of the targets, 
seems larger in the LMC. 

One puzzling fact becomes noticeable when the available Galactic samples are separated 
into luminous and low-luminosity YSOs. The YSOs in our LMC sample are luminous and 
thus should in principle be compared to the Galactic high-luminosity sample. However, it 
is apparent from Figures [HI and [9] that the properties of the CO2 profiles observed towards 
LMC YSOs more resemble those of low-luminosity Galactic YSOs rather than those of 
the high-luminosity Galactic YSOs, both in terms of water-rich fraction and concentration 
ratio. Admittedly the size of the LMC sample is small thus these results are tentative 
for now. We discuss how the temporal and spatial evolution of the physical conditions in 
molecular clouds can affect the observed properties of the different samples we compare. 

We also discuss the implications of our results in terms of the star formation conditions 
in molecular clouds in the LMC. The ISM of the LMC has lower metallicity than the 
Galactic ISM (Zlmc ~ 0.4 Z©). This affects the conditions in molecular clouds in the LMC 
in two parallel ways. On one hand there is less carbon and oxygen, and less dust. On 
the other hand, due to reduced shielding, the physical conditions in molecular clouds are 
harsher. 
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CO in the gas phase is less abundant, and the gasrdust ratio is larger by a factor 2 to 3 
in the LMC than in the Galaxy. Observations also suggest that CO cores are smaller in the 
LMC, the result of selective photodissociation of CO in the outer parts of GMCs. If water 
ice abundance is unchanged, these effects would conspire to reduce the ratio of CO2 ice to 
water ice in the LMC. This is clearly not what we observe. On the other hand, the twined 
facts that the UV radiation field is stronger and the average dust temperature is higher in 
the LMC could increase CO2 production, which would be consistent with what we see in 
our sample. 

Our results are tantalizing but not yet conclusive. We suggest that dust temperature 
could be responsible for observed differences between the LMC and the Galaxy. CO ice 
is more volatile than CO2 and the CO profiles can be used as an indicator of thermal 
processing if different components correspond to environments of different volatility. 
The suggestion is that CO ice is more sensitive to the relatively small dust temperature 
differences mentioned previously. Furthermore, the possible increases in UV photolysis and 
thermal processing are related to metallicity. We have obtained Spitzer-IKS observations in 
the SMC and complementary ground-based data in both LMC and SMC, with the goal to 
incorporate CO ice profiles in our analysis, and to extend these studies to the even more 
metal-poor environment of the SMC. 
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Table 1. Sample of embedded LMC YSOs, with groundbased opti cal and near-I i l pho t ometry. ?7-B\^/ a nd JHK 
fluxes are from the MCPS and IRSF catalogs, respectively (iZaritsky et al.ll2004l : iKato et al.ll2007l ). 



ID RA DEC F!7 Fb Fv Fj Fj Fh Fk 

(d ™ (^Jy) (^Jy) (^Jy) (^Jy) (mjy) (mjy) (mjy) 



IRAS 04514-6931" 






04 


51 


11 


47 


-69 


26 


46 


9 


58.6 ± 


10.0 


179.4 


± 


18 


7 


324.4 


±40.7 


194. 


3 


± 


19.5 





168 


± 





006 





396 


± 





018 





730 


± 





067 


SSTISAGBIC J045228 


65- 


-685451.3 


04 


52 


28 


65 


-68 


54 


51 


3 


5.0 ± 


1.8 


26.4 


± 


2 


4 


61.5 


± 


4.6 


Ill 


1 


± 


7.7 





068 


± 





007 





275 


± 





020 





748 


± 





041 


IRASF04532-6709 






04 


53 


11 


03 


-67 


03 


56 





34.2 ± 


3.7 


56.6 


± 


4 


1 


107.6 


± 


22.1 













112 


± 





009 





190 


± 





018 





323 


± 





027 


SSTISAGBIC J050354 


56- 


-671848.5 


05 


03 


54 


56 


-67 


18 


47 


7 


50.1 ± 


10.1 


55.9 


± 


7 


9 


53.7 


± 


4.5 


61 





± 


5.2 





475 


± 





009 


2 


034 


± 





019 


6 


658 


± 





061 


SSTISAGBIC J051347 


73- 


-693505.1 


05 


13 


47 


79 


-69 


35 


05 


1 



















































389 


± 





037 


SSTISAGBIC J051449 


41- 


-671221.5 


05 


14 


49 


41 


-67 


12 


21 


5 


5.7± 


2.5 


14.0 


± 


1 


5 


19.7 


± 


4.7 


73 


8 


± 


4.0 





224 


± 





008 


1 


112 


± 





021 


4 


811 


± 





133 


IRAS 05240-6809 






05 


23 


51 


14 


-68 


07 


12 


4 






68.2 


± 


10 


4 


780.4 


± 


108.5 










5 


936 


± 





711 


6 


021 


± 





388 


13 


158 


± 





727 


IRAS 05246-7137 






05 


23 


53 


95 


-71 


34 


44 


































055 


± 





004 





484 


± 





018 


3 


966 


± 





037 


MSXLMC464 






05 


24 


13 


30 


-68 


29 


59 





143.6 ± 


12.0 


192.8 


± 


18 


3 


255.1 


± 


23.7 










1 


765 


± 





032 


6 


909 


± 





127 


20 


231 


± 





373 


SSTISAGBIC J052546 


49- 


-661411.3 


05 


25 


46 


52 


-66 


14 


11 


3 


6.0 ± 


1.4 


15.5 


± 


1 


7 


34.2 


± 


2.0 


58 


8 


± 


6.6 





110 


± 





007 





462 


± 





017 


3 


254 


± 





060 


IRAS 05328-6827" 






05 


32 


38 


59 


-68 


25 


22 


2 































191 


± 





007 





028 


± 





418 


9 


111 


± 





168 


MSX LMC 1786" 






05 


37 


28 


07 


-69 


08 


48 









18.4 


± 


2 


4 


27.6 


± 


3.7 













254 


± 





012 





611 


± 





034 


1 


180 


± 





076 


SSTISAGBIC J054059 


29- 


-704402.6 


05 


40 


59 


29 


-70 


44 


02 


8 






31.1 


± 


6 


4 


26.6 


± 


3.8 


31 


9 


± 


3.5 















289 


± 





032 





780 


± 





044 


IRAS05421-7116'' 






05 


41 


25 


OS 


-71 


15 


32 


7 


50.5 ± 


4.8 


120. f 


± 


5 




57.4 


± 


7.8 













423 


± 





026 





690 


± 





044 


1 


285 


± 





058 


IRAS 05452-6924 






05 


44 


50 


23 


-69 


23 


04 


7 


12.9 ± 


2.5 


23.5 


± 


1 


5 


35.0 


± 


2.3 


78 


5 


± 


4.2 





181 


± 





007 





321 


± 





012 





671 


± 





025 



^sources from the Spitzer archive; all other sources are from the SAGE-Spec program. 
^JHK fluxes from the 2MASS catalog Jcutri et aLlbooj^ . 



Table 2. Sample of e mbedded LMC YSOs, with Spitzer IRAC and MIPS photometry from the SAGE catalog 
(IMeixner et al.l 120061 ) . The last column provides alternative identifications and known classifications from the 

literature. 



ID Fs.e,,,,, F4_5 j,„, F5.8,im Fg jim ^24^01 F70 Figo^im Comments 

(mJy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 



IRAS 04514-6931" 






5.72 ±0.494 


18 


90 


± 





547 


52 


70 


± 


1 


529 


96 





± 


6 


386 


2538.0 ±16 


590 


2538 





± 16 


590 


15820 ± 127.40 


1257 ±336.70 


MSX LMC 1200 


SSTISAGEIC J045228 


65 


-685451.3 


3.21 ±0.173 


6 


05 


± 





169 


10 


80 


± 





288 


23 





± 





514 


198 


3± 


1 


364 


217 


4 


± 


1 


064 


1435 ± 


24.24 


2793 ± 81.96 


YSO 


(1) 




IRAS F04532-6709 






2.03 ±0.350 


4 


00 


± 





171 


9 


10 


± 





381 


16 


9 


± 





880 


391 


3± 


2 


158 


350 


1 


± 


2 


214 


2454 ± 


29.70 


3815 ± 110.10 








SSTISAGEIC J050354 


56 


-671848.9 


24.50 ± 0.475 


35 


60 


± 





514 


45 


40 


± 





869 


64 


5 


± 


1 


579 


180 


6± 


4 


059 


307 


2 


± 


2 


703 


7374 ± 


72.46 


1091 ±385.10 


YSO 


(1) 




SSTISAGEIC J051347 


73 


-693505.1 




3 


13 


± 





158 


8 


25 


± 





417 


17 


6 


± 





830 


372 


5± 


1 


899 


276 


4 


± 


2 


949 








YSO 


(1) 




SSTISAGEIC J051449 


41 


-671221.5 


19.20±0.586 


34 


30 


± 





651 


51 


20 


± 





939 


71 


3 


± 


1 


238 


281 


7± 


1 


634 


294 


8 


± 


1 


674 


1243 ± 


19.92 




YSO 


(1), 


ST4 (2) 


IRAS 05240-6809 






22.50 ±0.685 


33 


40 


± 





529 


42 


20 


± 





860 


56 


8 


± 





797 


1216 


0± 


6 


113 


1270 





± 


7 


630 


3962 ± 


42.34 


1384 ± 91.71 


YSO 


(1), 


ST7 (2) 


IRAS 05246-7137 






34.50 ±0.694 


62 


50 


± 


1 


799 


92 


90 


± 


1 


672 


119 





± 


1 


877 


342 


6± 


1 


588 


348 


2 


± 


2 


042 


1490 ± 


19.09 




YSO 


(1), 


MIPS-SED |(3) 


MSX LMC 464 






63.60 ±1.890 


89 


50 


± 


2 


576 


122 


00 


± 


2 


684 


175 





± 


4 


752 


868 


0± 


4 


333 


668 





± 


4 


412 


6741 ± 


58.81 


1961 ± 421.60 


Hll? 


(4) 




SSTISAGEIC J052546 


49 


-661411.3 


30.00 ±0.668 


56 


60 


± 


1 


554 


79 


80 


± 


1 


430 


115 





± 


1 


958 


504 


7± 


3 


381 


498 


4 


± 


3 


021 


1417 ± 


34.56 


17510 ±363.90 


YSO 


(1), 


ST3 (2) ^ 


IRAS 05328-6827" 






77.24 ±3.294 


148 


20 


± 


3 


771 


212 


90 


± 


3 


233 


261 


7 


± 


7 


766 


1070 


0± 


5 


392 


919 


3 


± 


5 


637 


3052 ± 


63.57 


1113 ±230.80 


YSO 


(5), 


MIPS-SED '(3) 


MSX LMC 1786" 






11.20 ±0.500 


34 


20 


± 





784 


88 


10 


± 


1 


782 


142 





± 


2 


200 


1444 


0± 


9 


927 


906 





± 11 


460 


4139 ±231.90 


2002 ±636.10 


30Dor-17, MIPS-SED (3) 


SSTISAGEIC J054059 


29 


-704402.6 


4.51 ±0.202 


10 


30 


± 





297 


17 


20 


± 





480 


29 


5 


± 





581 


209 


7± 


1 


055 


216 





± 


1 


394 


945 ± 


12.78 




YSO (1) 




IRAS 05421-7116 






3.70 ±0.272 


4 


02 


± 





120 


7 


55 


± 





311 


13 


3 


± 





949 


331 


4± 


1 


770 


329 





± 


2 


454 


2040 ± 


22.66 




YSO (1) 




IRAS 05452-6924 






2.36 ±0.081 


3 


62 


± 





114 


8 


58 


± 





176 


19 


4 


± 





406 


674 


9± 


3 


158 


659 


3 


± 


3 


609 


2967 ± 


47.63 




YSO (1) 





^sources from the Spitzer archive; all other sources arc from the SAGE-Spcc program. ^ two epochs of 24-/j,m photometry arc available. 



References. — ("!') Iwhitnev et al.l J2008l t: (21 IShimonishi et al.l j200a t: f3l lvan Loon et al.l ^20091^ : M ^Kastner et allhOOat) : (5) \van Loon et al.M2005; 
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IRAS 05328-6827 



MSX LMC 464 




wovelength (;^m) 
IRAS 05246-7137 



wovelength (^m) 

MSX LMC 1786 




wovelength (;im) 
SAGEIC J052546.49-66141 1.3 



wovelength (fim) 

IRAS 04514-6931 




15 20 25 

wovelength {/im) 



15 20 25 

wovelength (/ioi) 



Fig. 1. — Spitzer-IRS spectra of the 15 sources discussed in this work. From left to right 
and top to bottom, sources are shown in decreasing order of their measured IRAC flux at 
5.8 //m. The insets in each graph show the spectral regions where ice bands can be found, 
namely the CO2 ice feature at 15.2 /im (lower right corner) and the 5 — 7/xm complex that 
includes the water ice O— H bending mode at 6/xm (upper left corner). 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 2. — Fits to YSO spectra using PAHFIT: spectra with all fitted components (upper 
panels) and continuum-subtracted spectra (lower panels). The total fit is shown (red) as 
well as the individual components of the fit: dust continuum (blue), PAH emission (green) 
and atomic and molecular emission lines (magenta) — in some cases there is a very small 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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6 8 10 12 14 16 18 20 5 8 10 12 14 16 18 20 

wovelength (^tm) wovelength (/^m) 



Fig. 3. — Continuum fits to YSO spectra without PAH emission. For each object the total 
continuum fit is shown (red dashed hne, low-degree polynomial) as well as the optical depth 
spectra. The two objects at the top have strong silicate absorption and weak CO2 ice; the 
two objects at the bottom exhibit either a self-absorbed or emission silicate feature and 
strong CO2 ice. 
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Table 3. Laboratory spectra use d in the COo ice fit s from 

and Wfiite et al. ( 20091 . 



Efirenfreund et al 



fll997l . E97) 



fl2009l . W09). 



mixture 


temperature 


reference 




(K) 




H2 0-rich polar ices 


H2U:CU2 = 100:14 


10 


E97 


H20:C02:CO = 100:20:3 


20 


E97 


CO- or CO2 


-rich apolar ices 




C0:C02 = 100:4 


10 


E97 


C0:C02 = 100:8 


10 


E97 


C0:C02 = 100:16 


10 


E97 


C0:C02 = 100:21 


10 


E97 


C0:C02 = 100:23 


10 


E97 


C0:C02 = 100:26 


10 


E97 


C0:C02 = 100:70 


10 


E97 


H20:CO:C02 = 1:50:56° 


10 


E97 


CO:02:C02 = 100:50:4 


10 


E97 


CO:02:C02 = 100:50:8 


10 


E97 


CO:02:C02 = 100:50:16 


10 


E97 


CO:02:C02 = 100:50:21 


10 


E97 


CO:02:C02 = 100:50:32 


10 


E97 


pure CO2 


10 


E97 


pure CO2 


50 


E97 


pure CO2 


80 


E97 


annealed ices 


H20:CH30H:C02 


5, 10,..., 135* 


W09 



''This mixture contains minimal traces of water. 
''5-K temperature steps. 
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polar HjO-rich ices 



^ 0.4 



o 0.6 




□ polor CO-rich ices 




15.0 15.5 16.0 16.5 14.5 15.0 15.5 16.0 16.5 

wavelength {fj.m) 

annealed ices 



0.0 



0.5 



1.0 



1.5 



2.0 



wavelength (junn) 
apolar pure CO2 ices 



1 1 




c6, c 


t iokI 


1 \ 


1 A ^ 

1 r 


CO2 c 


i 50kI 


~iv\ \ 
1 \1 V 
1 1 \\ 


I A] z'/ ^ 
' 1 r 






' \ \ ' 

1 \ V w 




CO2 c 


s bok' 


1 N\ ft 


1 










1 





14.5 15.0 15.5 16.0 16.5 
wavelength (jum) 




14.5 15.0 15.5 16.0 16.5 
wavelength (yu.m) 



Fig. 4. — Representative laboratory ice profiles for the different types of mixtures considered 
here: H20-rich, CO-rich, pure-C02 and annealed ices. The profiles used are listed in Ta- 
bleEl For the CO-rich ices the following compositions are shown: CO:02:C02 = 100:50:16, 
C0:C02 = 100:16 and C0:C02 = 100:70. Dashed (red) lines are the laboratory profiles 
smoothed to match the resolution of the Spitzer IRS LL mode (see text). The double- 
peaked structure of the pure-C02 and annealed ices is smoothed out to such an extent that 
it is no longer clearly identifiable. The vertical lines delimit the red and blue wing of the 
profile (see text). 
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14.0 14.5 15.0 15.5 15.0 16.5 17.0 14.0 14.5 15.0 15.5 16.0 16.5 17.0 

wavelength (/im) wavelength {f^m) 




14.0 14.5 15.0 15.5 16.0 16.5 17.0 14.0 14.5 15.0 15.5 16.0 16.5 17.0 

wavelength {f^m) wovelength (jum) 

-0.1 1 1 -0.11 J 

SAGEIC J052546. 49-66141 1.3 : 




0-5 F i 0-5 \r : 

14.0 14.5 15.0 15.5 16.0 16.5 17.0 14.0 14.5 15.0 15.5 16.0 16.5 17.0 

wavelength (/im) wavelength (/xm) 

Fig. 5. — CO2 profiles and best fits to the observed profiles (magenta), using a combination 
of either polar and apolar ices (left) or polar and annealed ices (right) — see text for ex- 
planation. The different components are color-coded as follows: water-rich (red), CO-rich 
(green), pure CO2 (blue) and annealed (cyan). Reduced values are also given. Objects 
are plotted in order of decreasing total column density from top to bottom (TableH]). 
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Fig. 5. — Continued. 
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Fig. 



5. 



, — Continued. 
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Fig. 



5. 



, — Continued. 
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Fig. 6. — CO2 features for the three YSOs for which no fits with laboratory spectra were 
performed. 



Table 4. Column densities and model fit information for the observed 15.2-/im CO2 features. Uncertainties are 
statistical and do not include those associated with the continuum determination. Fit components are water-rich 
(polar), CO-rich (apolarl), pure-C02 (apolar2) and annealed ices. For each component percentages (in terms of 
column density) are given. For the polar + apolar ice fits the best-fit C02:C0 ratio is provided; for the 
polar + annealed fits the best-fit annealed temperature is given, values for the fits are listed. Additional ice column 
densities are provided in the last co lumn for the 4.3-/im feature of CO2 and 3-/xm feature of water ice, either from th e 



literature (IShimonishi et al. 



20081 ) or r ecomp uted using published spectra (IShimonishi et al. 
2005al ) — see text and table notes. 



2008 



van Loon et al. 













P 


olar + apolar 


fit 








polar - 


- annealed fit 












ID 




N{C02)f5M™ 
10^" cm^^ cr 


polar 
% 


apolarl 

% 


apolar2 

% 


C02:C0 




polar 

% 


annealed 

% 


T 
K 


T range 
K 




best fit 


N(C02) 


4.3)jm 

10 1^ 


N<B20)3^,„ 
-JO 

cm 


IRAS 05240-6809 




12.69 ±0.28 


44.8 


60 ± 2 


20 ± 2 


19 ± 2 


> 1.12 


1.18 


54 ± 


1 


46 ± 1 


115 


110 - 120 


3.35 


apolar 


rr , + 4U 




~40'' 


IRAS F04532-6709 




9.91 ±0.46 


21.4 


66 ± 2 


34 ± 2 




> 1.12 


3.26 


66 ± 


3 


34 ± 3 


120 


115 - 130 


3.71 


apolar 








SAGEIC J052546.49 


-661411.3 


8.83±0.18 


49.9 


68 ± 4 


8± 4 


24 ± 3 


> 1.12 


3.74 


12 ± 


4 


88 ± 4 


110 


110 - 115 


3.09 


annealed 


15 5+^3- 


i 


29 7+* '> 


SAGEIC J051347.73 


-693505.1 


8.37 ±0.93 


9.0 


59 ±14 


18 ±14 


23 ±13 


> 1.12 


0.64 


42 ± 


7 


58 ± 7 


115 


100 - 130 


0.53 










IRAS 04514-6931 




7.99 ±0.08 


98.3 


78 ± 3 


9± 3 


13 ± 2 


0.70 ±0.47 


10.01 


73 ± 


1 


27 ± 1 


125 


110 - 130 


8.21 


annealed 








MSXLMC 1786 




7.82 ±0.47 


16.6 


54 ± 3 


27 ± 4 


19 ± 3 


> 1.12 


3.90 


40 ± 


2 


60 ± 2 


115 


110 - 125 


3.28 


annealed 








SAGEIC J051449.41 


-671221.5 


6.84 ±0.47 


14.7 


51 ± 6 


44 ± 6 


5± 8 


> 1.12 


2.08 


49 ± 


3 


51 ± 3 


120 


110 - 130 


2.29 




„ , + 4.3 






IRAS 05421-7116 




4.95 ±0.86 


5.8 


69 ± 6 


31 ± 6 




> 1.12 


2.92 


60 ± 


8 


40 ± 8 


115 


110 - 130 


2.74 










SAGEIC J045228.65 


-685151.3 


2.65 ±0.44 


6.0 




73 ± 22 


27 ±22 


0.70±0.57 


3.72 






100 


125 


115 - 130 


3.61 










IRAS 05452-6924 




2.36 ±0.42 


5.7 


55 ± 14 


20 ±20 


25 ±10 


0.10±0.29 


0.62 


29 ± 


9 


71 ± 9 


135 


135 - 135 


0.58 










IRAS 05328-6827 




0.90±0.15 


5.9 


55 ± 8 




45 ± 8 




0.74 






100 


110 


110 - 110 


0.83 








^■^°~0.27 


SAGEIC J050354. 56 


-671848.5 


0.87 ±0.22 


4.0 


50 ± 8 




50 ± 8 




1.47 


46 ±10 


54 ± 10 


120 


110 - 130 


1.44 










SAGEIC J054059.29 


-704402.6 


0.63 ±0.77 


0.8 
































IRAS 05246-7137 




0.46 ±0.18 


2.5 
































MSX LMC 464 




0.45±0.17 


2.7 

































assuming A — 1.1 X 10 cm molecule "'" jcerakines et aljllQQs! ) and integrating between 14.4 and 16.8 ^im. 

estimated from Figure 1 by lshimonishi et al.| J2008| ) . new estimate using the spectrum published by Ivan Loon et al.| J2005a| ) . 



00 



SAGEIC J052546-66141 1.3 



5.5 



6.0 



6.5 



7.0 



7.5 




5.5 



6.0 



6.5 

wavelength (yum) 



7.0 



7.5 



Fig. 7. — Optical depth in the 5 — 7 fim region, respectively for 

SAGEIC J052546.49-661411.3 and IRAS 05240-6809. The filled circles indicate the 
position of unresolved H2 lines, of which two are clearly present in the spectrum 



of IRAS 05240 



( jHudgins et al. 



680 9. We also show typical laboratory ice models for pure water 



19931 ). scaled to match the water column densities measured using the 3-/im 



feature (TableHl), taking into account their error bars (shaded areas) . The dashed lin e is an 



Zasowski et al. 



(120091 ) 



optical-depth spectrum for the Galactic YSO IRAS 04016+2610 from 
showing water ice and the additional ice feature at 6.8 fim (spectrum scaled by a factor 1.5). 
Both LMC objects clearly show the presence of water ice and SAGEIC J052546. 49— 661411.3 
also shows the 6.8-/xm feature (see text). 
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Fig. 8. — Fraction of CO2 locked in the water-rich (polar) component, for the LMC sam- 
ple (top) and Galactic YSOs and background so urces (bottom). Galactic d ata come from 



Gerakines et al. 
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Fig. 9. — CO2 ice versus H2O ice col umn densities. W e comp iled mea surements from the liter- 



ature for Galactic sou r ces (c ircles) 



Gerakines et al 



pink). 



Zasowski et al. 



(I2OO7I . 



200% and 



2009 



Knez et al. 



cyan 



1999 



Pontoppidan et al. 



J gre en) 



Nummelin et al 



(12008 



, gold), and 



(2001 



Whittet et al. 



(120051 ) (black). Filled circles are YSOs while open circles 



are background sources (i.e. quiescent s ightlines). Star symbols are the seven LMC mea- 



surements: from iShimonishi et al 



(12008 



blue stars) and our CO2 ice measurements com- 
bined with literature H2O ice measurements (red stars). The full lines are median ratios 
N(C02)/N(H20) for Galactic high-luminosity and background sources (~0.2, green), Galac- 
tic low-luminosity sources (~0.3, gold) and LMC sources (~0.32, red) — see text for full 
explanation. 



